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■ We solve the magnetohydrodynamic (MHD) equations governing axisymmetric flows 

around neutron stars and black holes and found all possible solution topologies for adi- 
abatic accretion. We divide the parameter space spanned by the conserved energy and 
angular momentum of the flow in terms of the flow topologies. We also study the possi- 
bility of the formation of the MHD shock waves. 
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^ ' 1. Introduction 

^ ' In recent years, the study of magnetohydrodynamic (MHD) accretion flow around 

compact objects has become very important since the magnetic field in ubiquitous 
in the universe and it should play a part in the black hole accretion, especially 
^ ■ close to the inner disk. The pioneering works of Mestel (1967)^ and Weber & Davis 

^ ! (1967)2 etc. were further extended by Chakrabarti (1990)^ for rotating compact 

I stars and obtained a few global solutions for accretion and winds. The study of 

■ MHD flow in Kerr geometry is also carried out by Takahashi et al. (1990) .'^ They 

discussed the properties of smooth trans- Alfvenic flows and focused on the possible 
' ways to extract energies out of the BH. Nitta et al. (1991)^ presented an analytical 

, solution of general relativistic Grad-Safranov equation around a rotating BH and 

obtained few solution topologies. So far, a fully self-consistent study of the global 
^ ' solutions around a compact star in terms of flow parameters was not explored. In this 

, paper, we wish to provide a complete study of the trans-magnetosonic accretion flow 

5h ' around compact objects in presence of both the radial and toroidal magnetic fields 



using Paczyhski-Wiita^ pseudo-Newtonian potential. We discuss how the nature 
of the solution topologies depend on the parameter space spanned by the energy 
{£) and the angular momentum (L) of the flow. We identify eighteen types of 
solution topologies and divide the parameter space accordingly. We also flnd that 
the standing magneto-hydrodynamic shock wave forms which could give rise to 
the high energy particles ( few MeV) through the shock acceleration mechanism of 
the electrons. These accelerated particles produce power-law synchrotron radiations 
which could explain some of the spectral properties of BH candidates.^ Moreover, 
soft photons are energized by the hot electrons of post-shock flow through inverse 
Comptonization. In this paper, we provide all possible solutions of MHD flow with 
or without shocks and classify the parameter space according to the nature of the 
flow solution. 
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2. Behaviour of MHD shocks, flow solutions and parameter space 

An example of complete MHD solution with standing, slow-magnetosonic shock is 
presented in Fig. la. We consider the Alfven velocity and the radius to be i?o = 10^° 
cm s~^ and Va = let'em respectively. The central mass was chosen to be IOMq 

for illustration purpose. The flow parameters arc (£, L) = (1.5, 1.45). The shock 
location is predicted at = 1.153 and denoted by the vertical dashed line. The 
solution towards the BH is depicted by a single- arrowed curve while the sub-Alfvenic 
double-arrowed c;urve is appropriate for a shock free neutron star solution. The 
dash-dotted vertical line with triple-arrow represents a typical shock transition at 
the boundary of the neutron star. The hot and dense post-shock flow is the most 
important region since it is responsible for generating the hard X rays. 

In Fig. lb, wc draw the parameter space in £-L plane and divide it into various 
regions according to the nature of the solution topologies. We obtain the curve ADE 
consideriner the snecial case where the radial velocitv at all the magnetosonic points 




(c) (d) 

Fig. 1. (a) A complete solution topology with standing magnetohydrodynamic shock wave.* (b) 
Division of the parameter space according to the nature of solution topologies.* (c-d) MHD flow 
solutions for different (8, L) pair. Radial distance is plotted along horizontal axis while radial 
velocity is plotted along vertical axis. Greek alphabets mark various types of solutions drawn with 
parameters from different regions marked Fig. lb.* 
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vanishes. The region shaded with dot-dashed Unes is forbidden for the flow solution. 
The curve PG is obtained when flow velocity is identical to the Alfven velocity at 
the magnetosonic point. The regions bounded by PBEGHP, ABCFGUA, RDEQR 

and QGUQ possess Bondi-likc slow, Rotational-slow, Rotational-fast and Bondi- 
like fast magnetosonic points respectively. The parameters from the region shaded 
with dots (SKTS) allow magnetosonic shock transition. The region surrounded by 
PCKOTUQP may produce non-steady shock sincx; shock conditions are not satisfied 
but entropy at the inner sonic point is higher than the outer sonic point. 

We further sub-divide the parameter space according to the solution topologies. 
We identified eighteen distinct types of solutions as shown in Fig. l(c-d). The ac- 
cretion solutions for BH and neutron star differs only through the inner boundary 
conditions. Therefore, for a set of flow parameters, one needs to assign the proper 
boundary conditions to decide whether the solution is for BH or for neutron star. 
We marked each solution by a Greek alphabet a, (3,.. etc which denote the pa- 
rameters in £-L plane (Fig. lb) for which the solutions are drawn. In Fig. l(c-d), 
the Bondi-like (slow/fast). Rotational slow and Rotational fast magnetosonic points 
are denoted by mi, m2, and ms respectively. For L = 0, the solution reduce to the 
Bondi flow solution. 

3. Concluding Remarks 

We identified all possible solutions including those containing standing shocks in 
a MHD accretion flow around compact objects and divided the parameter space 
according to the nature of flow solutions. We identified another region where the 
flow may exhibit non-steady shocks very similar to those in a hydrodynamic flow. 
The post-shock flow inverse Comptonizes soft photons either from the cooled Kep- 
lerian disk or from the synchrotron radiation and re-emits them as hard X-rays.^ In 
addition, high frequency QPOs of hard X-ray could be obtained since MHD shocks 
form closer to the BH horizon. 
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